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Abstract 

The idea of defining a gravitational singularity as an obstruction 
to the dynamical evolution of a test field (described by a PDE) rather 
than the dynamical evolution of a particle (described by a geodesics) 
is explored. In particular, the concept of field regularity is introduced 
which serves to describe the well-posedness of the local initial value 
problem for a given field. In particular this is applied to (classical) 
scalar fields in the class of curve integrable spacetimes to show that 
the classical singularities do not interrupt the well-posedness of the 
wave equation. 


1 Introduction 

One of the biggest surprises that General Relativity (GR) has given us is that 
under certain circumstances the theory predicts its own limitations. There 
are two physical situations where we expect General Relativity to break 
down. The first is the gravitational collapse of certain massive stars when 
their nuclear fuel is spent. The second one is the distant past of the universe 
when the density and temperature were extreme. In both cases we expect 
that the geometry of spacetime will show some pathological behaviour. 

The first step towards a mathematical characterisation of the circumstances 
under which GR breaks down was achieved in the seminal work of Penrose 
and Hawking in their singularity theorems (see [I], Chapter 8). The general 
structure of the theorems establish that if on a spacetime ( M,g a b ): 

• a condition on the curvature 

• an appropriate initial or boundary condition 
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• and a global causal condition 

are satisfied then ( M,g a b ) must be geodesically incomplete |2]. 

The characterisation of singularities in terms of geodesic incompleteness 
requires us to consider spacetimes with metrics of differentiability at least 
C 1,1 (also denoted C 2 ~) in order to have a well-defined notion of geodesics. 
For any metric with lower differentiability one does not have uniqueness of 
the geodesic equation so it is not clear how one can define a precise notion 
of singular behaviour using point particles as probes to test the geometry. 
This was the motivation for Clarke to introduce the notion of a □- globally 
hyperbolic spacetime [3] which involves probing the geometry of spacetime 
with classical fields rather than point particles. 

Earlier work by Wald [3] gave a prescription to define dynamics in static 
non-globally hyperbolic spacetimes. Based on similar techniques, Kay and 
Studer jSj determined the boundary conditions for quantum scalar fields on 
singular spacetimes with conical singularities representing cosmic strings. 
Subsequently Horowitz and Marolf [6j, used Wald’s approach to study the 
theory of quantum free particles in static spacetimes with timelike singu¬ 
larities. They used the term quantum regular if the evolution of any state 
is uniquely defined for all time. The main technique used for this was to 
notice that that if the the spatial Laplace-Beltrami operator is essentially 
self-adjoint in L 2 (E / , \/hd 3 x) (where T,' is a three dimensional geodesically 
incomplete manifold and \Jhd 3 x the volume form of the induced metric on 
S') then using standard properties of self-adjoint operators, a unique evo¬ 
lution of the wave function is obtained. Moreover, the classical singularities 
disappear in the sense that there is no freedom in the boundary conditions 
to define the state evolution. Later work, by Ishibashi and Hasoya [7] used 
similar techniques to investigate the evolution of the wave equation O g 4> = 0 
in static singular spacetimes by focusing on changing the function space from 
L 2 (E / , y/hd 3 x) to iL 1 (E / , \/hd 3 x). The main reason for this is that finiteness 
of the energy states implies the finiteness of the H 1 norm. In addition, ex¬ 
ploring different Hilbert spaces is important for quantum theory in curved 
spacetimes. Differences in the corresponding quantum field theories might 
provide useful insight into understand the behaviour of quantum states near 
singularities in curved spacetimes. Finally, Ishibashi and Hasoya used the 
term wave regular if the initial value of the wave equation has unique solu¬ 
tions in the whole spacetime with no arbitrariness in the choice of boundary 
conditions. Vickers and Wilson |5] also studied the problem of conical sin¬ 
gularities from Clarke’s perspective and Wilson in [9j showed that one could 
obtain dynamic evolution subject to constrains on the initial data and a flux 
condition in the singularity. 

The link between the concept of □- globally hyperbolic, quantum regular 
and wave regular is to redefine a singularity in spacetime not as an obstruc¬ 
tion to geodesics or curves but as an obstruction to the dynamics of test 
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fields. Nevertheless, each concept has it own characteristics. While quan¬ 
tum regularity probes spacetime with a quantum free particle, the notion 
of □- globally hyperbolic and wave regular uses the classical wave equation. 
Also, while quantum regularity and wave regularity look to singularities in 
terms of boundary conditions, a □-globally hyperbolic approach identifies 
the singularity as an interior point in a spacetime with low differentiability. 
In an heuristic manner, one can refer to field singularities as any approach 
to identify and characterise gravitational singularities as an obstruction to 
the evolution of test fields. This is in contrast to the standard approach 
where one uses geodesic incompleteness (which describes an obstruction to 
the evolution of a test particle) to identify singularities 
In this paper the particular case we will deal with is that of the wave equa¬ 
tion and we define a point p in M to be strongly wave regular if there is a 
neighbourhood of p in the spacetime (M, g a b) such that there is a lens-shaped 
domain U containing p (see figure 1) and there is a triple ( V , Q, 1Z x S) such 
that the initial value problem for the wave equation \H g <f> = f on U + is 
locally well posed in the following sense: 

• There exists a solution in the function space V(U + , y/—gdfix). 

• The solution is unique in the function space Q{U + ,^/ zr gd 4: x). 

• The solutions in the space Q(U + , ^fi^^gd^x) depends continuously with 
respect initial data in function space 7H(Eo, Vhdfix) x 5(Xo, Vhdfix). 

A weakly wave regular point p only satisfies the first two conditions. A 
strongly wave regular spacetime is defined to be one such that every point p 
in (JA,g a b) is strongly wave regular. 

In this paper we refine and provide full details of Clarke’s arguments in [31 
where he outlined this issue for what he called curve integrable spacetimes 
and give a detailed description of the techniques needed to show that solu¬ 
tions to the wave equations exist and are unique. Roughly speaking a curve 
integrable spacetime is one in which the integrals of both the connection 
and curvature along a curve are bounded (see condition (4) of the main the¬ 
orem for a precise description). From a physical point of view a spacetime 
is a curve integrable spacetime if there is a set C that defines a range of 
timelike directions which are transverse to any shock or caustic that may be 
present. The theorem proved by Clarke in [31 required both the quadratic 
and linear part of the Riemann tensor (in terms of the Christoffel symbols) 
to be separately integrable along the timelike directions. However Clarke 
also conjectured that one can prove existence and uniqueness of the wave 
equation if one only required the weaker condition of the integrability of the 
Christoffel symbols and the Riemann tensor. In the present paper we prove 
this conjecture. In addition, we also show continuity with respect initial 
data and establish strong wave regularity. 
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2 Curve Integrable Spacetimes 


In this section we establish the well-posedness of solutions to the initial 
value problem for O g 4> = f where D 9 is the wave operator given by a metric 
g a b that corresponds to a curve integrable spacetime. The proof is in four 
main steps. The first step is to define an energy inequality which is an 
inequality between an integral of the function and its derivatives at a future 
time bounded above by an integral of the function and its derivatives at 
the initial time and the source function. The second step is to show self¬ 
adjointness of the wave operator in some appropriate function space and to 
use the Hahn-Banach theorem to prove existence of a solution. The third 
step is to notice that the energy inequality allows us to show uniqueness 
of our solution. The last step is to conclude that the solution depends 
continuously on the initial data again using the energy inequality. 

Following Clarke J3j, we introduce an enlarged notion of a solution for these 
spacetimes. The general geometric background we use to define a generalised 
solution is the existence of a lens-shaped domain U properly contained in a 
open subset P of a 4-dimensional spacetime (M,g a b). 

This means that there is a smooth map © : T x (—a, a) —> M where Sc M 
is a compact, C 1 co-dimension 1 sub-manifold with boundary, with the 
property that: 

• ©(•, 0) : T —>• So is the identity map. 

• 0(x, r) = 0(x, s) for any x E dT, t,s E (a, b). 

• for any fixed s E (a, b), 0(5, s) is an 3- dimensional spacelike hyper¬ 
surface. 

• away from dT x (a, b), 0 is a diffeomorphism. 


We denote the region from 0 to r > 0 by U '+ and from t' < 0 to 0 by U~,. 
Notice that given coordinates x l on T, 0 provides coordinates (t,x 1 ) for 
the region away from the image of dT,. We will therefore always choose 
charts such that the time coordinate coincides with the time coordinate 
given by 0 when working in coordinates. 
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Figure 1: The general geometric setting. 

We start by considering the case of a C 2 solution to the problem: 

= / ( 1 ) 

in the region U+ where g a b is a smooth metric with C 2 x C 1 initial data 

(¥h 7r ) 


0|s o = ¥> 

90, 

9f 


So = ^ 


( 2 ) 


(3) 


Then multiplying {[l]) by a test field p E (M+ ) and integrating by parts 

gives: 


= - j u} pt^x- j^4hd\ (4) 

0|s o = 9? (5) 

Clarke then generalised the idea of a solution to 0 ><© and ([3]) in the region 
U+ for a low differentiable metric to be a function cf) that satisfies condition 
0 V/9 E V(U+) with initial value 0. Notice that the expression is well 
defined for </> E H 1 (Z/+ , yf—gd^x) , </? E H 1 (T,Q,Vhd 3 x), ir E L 2 (T,Q,Vhd 3 x), 
f E L 2 (U+, ypgd^x) and <7 U E L 2 {U+ , y/ zr gd 4 x) where ypgd 4 x is the 
volume element of the metric g a b and Vhd 3 x is the induced volume given by 
the induced metric h a b on So- To simplify notation we define u g := ypgd 4 x 
and ph '■= Vhd 3 x. 
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2.1 The main theorem 

The theorem we will prove can be stated precisely as: 

Theorem 1 Let (M,g) be a Lorentzian manifold and p a point in an open 
subset Ll C M with compact closure such that there is a lens-shaped domain 
U satisfying: 

1. The components gij and g~^ 1 are C° ; 

2. The components gij are C 1 in AL\J + (p) 

3. gij,k exist and are L v loc [M) with p > 2; 

4■ there exist functions Rj kl which, interpreted as distributions, coincide 
with the Riemann tensor defined distributionally from g and gij t k- 

5. there is a non-empty set, C Cl 4 , and positive functions M, N : M + —>• 
M + such that, if 7 is a curve with d'y/ds £ C for all s then 


7 is future timelike; 


the integrals 

ra 



I 7 (a) := j o |rj fc (7(s))|ds 

(6) 

and 

PQj 



M a ) ■= ) o \Rjkih(s))\ds 

( 7 ) 

(where T is defined using the weak derivatives of g) are 

conver- 

gent, with 

/ 7 (o) < M(a), J 7 (a) < N(a) 

(8) 

and M{a),N(a ) 

—> 0, as a — > 0. 



then p is strongly wave regular. 

2.2 Energy inequality 

We follow the classical method of proving well posedness of the wave equa¬ 
tion by using an energy inequality as shown for example in Hawking and 
Ellis [Tj , Clarke [3] and Wilson |9J. The energy inequality gives an integral 
of the function and its derivatives at a future time bounded above by an 
integral of the function and its derivatives at the initial time and an integral 
of the source function over the region between the initial time and the future 
time. We will first assume that g a b is smooth. Then at the end we will give 
the extra requirements that the metric must satisfy in order that the energy 
inequality is still valid when the differentiability is below C 2 . 
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Consider </> a solution to the wave equation that is C 2 with energy-momentum 
tensor T ab given by: 

T ab = (s-V" - ^“V") P) 

Now choosing a smooth timelike vector field, Y a , (which in the theorem will 
be chosen to be the 4-velocity of a family of observers) we define the energy 
integral: 


E(t) = f T ab T a n bf i h (10) 

where n a is a future pointing vector normal to £ r . 

Then we use the divergence theorem on the domain lit : 

/ div (T a6 T a ) v g =[ T ab T a n b ^ h (11) 

JU T J oU T 

The left hand side takes the explicit form: 

J u+ (/^ Ta ) [/ - $ + T ab V b T a NVhd A x ( 12 ) 

where f = \H g 4>, N is the lapse function which satisfies Ndt = n a and \/h is 
the induced metric in the hypersurface So- 
The right hand side then takes the form: 



Now we introduce the following norms: 



The main reason for introducing these norms which are naturally related to 
Sobolev type norms is the following result which can be found in [9j. 

CiE{T)<mk?<C2E(T) (17) 


for constants C\. C 2 > 0. 
Also notice that: 
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Ml# ^ k 


dt 


(18) 


where k is a constant that depends on g a b- 

Now we can obtain the following bounds for all the terms in (12): 


g v^r i )f N g h )dt<K l m\\uf + 


u t ?) ( 19 ) 


d<p 




/ o \JT, t \ dxi 
where K \, K 2 are constants that depend on g a b and Y a . 


( 20 ) 


The last term in (12) is bounded by: 


iui 


T ij VjTi < K 3 (\ml+) 


( 21 ) 


where K 3 is a constant that depend on g a b and Vf,Y a . 

Estimating all the terms in ( |11[ ) by the bounds available (19), (20) and (21) 
gives the inequality: 

E(t) < £( 0 ) + M||/||“ + ) 2 + ki(ll <+) 2 ( 22 ) 

where ko, k\ are positive constants that depend on the metric g a b, the vector 
field T a and the covariant derivative VfcT a . 

Now rewriting (22) using © and (|18[) as : 


E(t) < E( 0) + fe 0 (||/||° + ) 2 + k 2 J E(r')dr' (23) 

Using Gronwall’s inequality the desired energy inequality is obtained: 


E(r) < K 4 (E(0) + (||/||°+) 2 ) for all r < h (24) 

where K 4 positive constant that depends on the chosen finite time t \.the 
metric g a b, the vector field T a and the covariant derivative Vf,T a . 

In term of the Sobolev norms we obtain the expression: 

(WU 2 < ((Wto) 2 + (ll/lfid 2 ) (25) 


where we say a < b if a < kb for some constant k. 

We have proved all these results for C 2 functions but we can readily extend 
them to the subspace of solutions of (jlj) in E[ l {U+ ,u g ) with all derivatives 
in L 2 {Ti r ,gh) for all r E [0,ti]- 
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We now look for the conditions required to obtain again the energy inequal¬ 
ities in the low differentiable setting. The basic requirement we need is that 
the we can apply Stokes’ Theorem. To our knowledge the optimum results 
are given by the following theorem: 


Theorem 2 Let Ll be a compact set with compact closure with Lipschitz 
boundary and let Z a be a vector field on an (n + 1)- dimensional manifold 
M with continuous metric g a b and metric volume element v g . If 

Z a € w£(M), 9ab G W^ + \M), 

then the Stokes identity holds: 



for j = iz a Vg and where dj = div(Z a )u g 


A proof of this can be found in mu. 

Unfortunately, the result can not be applied directly to our case to obtain 
the energy inequalities as it would require us to take G H? oc (M). The 
way to proceed is to again first assume G C 2 and look for the conditions 
needed on the metric and the vector field T a such that we obtain again the 
energy inequalities. 

Assume now that the metric and its metric volume element are continuous. 
This is enough to satisfy the hypothesis and allows us to apply Stokes’ 
theorem. Of course we would like to have the same expression as in ©• 
This requires the existence of a metric connection, i.e. V g a b = 0. In m , it is 
stated that sufficient conditions for the existence of a Levi-Civita connection 
are the existence of a connection V G Lf oc (M, v g ) and that g ab G L^ C (M). 
We say V G Lg oc {M , v g ) if Vx Y a G Lf oc (M, v g ) for any pair X a . Y a of C°° 
vector fields. If the Christoffel symbols satisfy T l - k G L 2 oc (M, v g ) then they 
define a Levi-Civita connection. For example, the class of Geroch-Traschen 
class of metrics satisfy the above conditions. 

It can be seen by direct inspection that the other inequalities require only 
that g ab G Lff c (M ) and g a b G Lff c (M ) which are also enough to maintain 
the results © and ( |18[ ). Finally the conditions on the timelike vector field 
Y a which is needed is that the covariant derivative is essentially bounded. 
The sufficient analytical conditions which guarantee the existence of such 
vector field are a key part of our result and are established in Proposition [2] 
shown in the appendix. 

Now we have all the analytical conditions needed to recover the energy 
inequality ( |25[ ) we can again extend the result to the subspace of solutions 
of (fl]) in // ( Ufi , Vg) with all derivatives in L 2 (S r ,/r/i) for all r G [0,fi]. 

For clarity we state the result as a Lemma: 
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Lemma 1 If (M,g) is a Lorentzian manifold and there is subset fl with 
compact closure such that there is a lens-shaped domain, Uff, in fl satisfying: 

1. The components gij and g~^ are C° ; 

2. The components g tJ are C 1 in M\J + (p); 

3. weak derivatives gtj^ exist and are L 2 oc (M, v g ); 

4 ■ there exist functions Rj kl which, interpreted as distributions, coincide 
with the Riemann tensor defined distributionally from g and gij,k- 

5. there is a non-empty set, C Cl 4 , and positive functions M, N : M + —>■ 
M + such that, if 7 is a curve with d'y/ds £ C for all s then 

• 7 is future timelike; 

• the integrals 

^ 7 ( 0 ) : = J o \ r )k('y( s ))\ ds ( 26 ) 

and 

J^a) := J* \I? jkl ('y(s))\d8 (27) 

(where T is defined using the weak derivatives of g) are conver¬ 
gent, with 

I~/(a) < M(a), J 7 (a) < N(a) (28) 

and M(a),N(a) —> 0, as a—>■ 0. 

then 



for all f> £ with all derivatives in L 2 (S r , ph) for all r £ [0,fi] 

and [——1 g cf) = f . 

2.3 Self-adjointness and Existence 

The next step required to establish the theorem is the self-adjointness of the 
operator in an appropriate function space. 

Consider the ^(U^Vg) norm. 



We shall single out two important subspaces of H 1 , u g ): 
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Vs 0 = C°°(n) s. t. V'lso = n a if, a |So = 0 and U g if £ L 2 (U t +, v g )} 

Vs tl = {w G C 00 ^) s. t. o;| Sti = n a u:,a |s T = 0 and □gV’ G L 2 (U^, v g )} 

Then the condition that g VJ is C° and that Z a = £ C°°(£l) and Z a = 

ifuj a G C°°(Ll) is enough to apply Theorem [2] and obtain: 


/ g ab if b u a v g + 

J K 



LOUglfVg 



g ab if b u a v g + 



ifUgUVg 


0 

0 


(31) 

(32) 


so combining the above equations give 


and we obtain 




(33) 


{Dgi/j,uj) = {ip,D g uj} (34) 

which proves the self-adjointness of D s for if £ Vs 0 , oj £ Vs t . 

In the following lemma we establish a density result needed later . 


Lemma 2 The spaces Vs 0 , Vs t are dense in L 2 (Uff,v g ). 


Proof of Lemma. It is enough to prove that Cff (lAff ) is contained in both 
sets.pTSJ 

Let (f £ The only thing we need to prove is that Z\ g <f £ L 2 {U^, v g ). 

Explicitly in coordinates we have that 


Ogcf = g 




d<f) 


+ r A 

1 7. 


d<f) 


dx'dxJ lJ dx k 


That the first term is in L 2 (U^,v g ) follows directly from g^ £ C° and 
4> £ C?(U+). The second term is in , v g ) because 



< oo (36) 


where we have used the fact that the connection is in L 2 oc (M, v g ) and 
Holder’s inequality. This allows us to conclude then that O g (f £ L 2 (U ^) 
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hence 4 > is in Vv 0 , U St| for any £ C£° (U^f ). □ 

The proof of existence uses the Hahn-Banach Theorem and the Riesz Rep¬ 
resentation Theorem (see [T2j p 199). We start by proving that □ 5 Vs 0 is 
dense in L 2 (£7^, v g ). It is enough to prove that 

□sV'eo = {x E L 2 (U+,v g )\(x,y) for all y £ D g V^} = 0 . 

Suppose 77 £ , because Vz t is dense in L 2 {U^. u g ), there exist a 

sequence {w n } in Vjv converging to rj in L 2 (U^f. v g ). 

Now we have: 


(^□jW n ) = (n g <j>,u n ) ->■ {O g (j),ri} = 0 (37) 

for all <p £ Vs 0 - Since, Vs 0 is dense, this implies that (||n s a; n ||° + ) 2 —> 0. 

lA r 

In order to continue the proof we need the result of the following Lemma. 

Lemma 3 For any element £ in So we have the following form energy in¬ 
equality: 

(kill,) 2 S (IP a «l£+) 2 (38) 

for 0 < r < t\. 


Proof of the Lemma. We prove the energy inequality by choosing a future 
pointing normal vector field. However, we can also repeat the argument in 
an identical form by choosing a past-pointing normal vector field to 
In that case we obtain the energy inequality: 


(WU 2 <((Wk 1 ) 2 + (IP 9 <+) 2 ) (39) 

for 0 < r < t\. 

Then energy inequality takes the form: 

(HIU 2 <(IIQ ,< + ) 2 On) 

for 0 < r < t\ and u> £ V£ t] . □ 


Using now the energy inequality (39) with (18) and integrating both sides 
from 0 to t± we obtain: 


(iwiLks /“(ip^„ifi + ) 2 <iT (4i) 

1 Jo T 

Hence, (||u ; n || 1 + ) 2 —>• 0 which implies using the Sobolev embedding theorem 
that (||cuJ| 7 °,+ ) 2 —t 0 hence r] = 0 . 
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We define the functional kf(O g u) = (/, w). We show the functional is 
bounded so we can apply Riesz’s Theorem. In that way kf define an el¬ 
ement G L 2 (Uv g ) such that kf(JD g uj) = (\k, O g co) 

Now using the energy inequality 


40 and Cauchy-Schwartz we obtain: 


k(n g u) = (f,u>) (42) 

< (/,/)(w,w) (43) 

< (\\f\\° u+ ) 2 m g u\\° u+ ) 2 (44) 

The Hahn-Banach theorem together with the fact that DgVso is dense in 
L 2 (^+, v g ) allows us to extend the functional to the whole L 2 (U^, v g ) with¬ 
out increasing the norm.Hence, kf is bounded in L 2 (U ' t +, v g ). 

Then Riesz’s Theorem implies there is a G L 2 (Uv g ) such that 

/ ^ = , f uu g 

Ju+ Ju+ 

for all oj G Vs t| D Then ^ is a solution of □ 9 \I / = /. 

To add initial conditions we construct a specific function q that satisfies the 
required initial conditions. We show this in detail in the following proposi¬ 
tion: 


Proposition 1 


There is a function (j> in L 2 (Uv g ) that satisfies: 


•a g fi = f 

• &ko = <p 

m d<t>\ _ 

• gtko = * 

with ip, 7T G H 2 (Tj o, Hh)nW 1,2r (Eo, ph) where the connection is in L 2p (Uk, v g ) 
and j + l = 1 


Proof of Proposition 1. Let q = ip + tn G H 1 (Uf)) which implies q|s 0 = p 
and || | s 0 = 7r - Then 


O g q = gV 


dq 


+ r A 

1 7. 


dq 


dx l dxi lJ dx k 


In order to that expression be in L 2 (Uk, v g ) we need to look for sufficient 
analytic conditions. The first term imposes that 9 Jg x j G L 2 (U^,v g ) so 
q must be in H 2 (U^,v g ). To analyse the second term we use Holder’s 
Inequality as follows: 
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where - + - = 1. 

p y 

If the connection only satisfies being in L 2 (Z00, u g ) then we need that q E 
If on the other hand we know that the metric is II /1,oc (Z^) 
then the condition q E H 2 (U^ , u g ) can be maintained. Hence, sufficient 
conditions for the initial data are that <p, it are in i7 2 (£o, ///JnTU 1,2r (Xo, Hh )■ 
Then applying Hahn-Banach theorem to the functional &/_□ 9 we obtain a 
function Ik' such that 

f + *' n gU v g = [ (/ - 

Ju t 1 Ju t 1 

which again satisfies = ^ |s 0 = 0 

Then the desired solution is given by 0 = ik' + g.Q 

That there are solutions in H k (U^,Vg) require higher energy inequality es¬ 
timates. The main steps are to rewrite the wave equation as a first order 
system, obtain energy inequalities in Sobolev Spaces with negative integer 
and then apply standard dualities. (See 11 -ij . Chapter 28: m, Chapter 7). 


2.4 Uniqueness and continuity with respect initial data. 

The proof of uniqueness follows directly from (25). Take two functions 0i 


and <f >2 in 7/ 1 , u g ) such that solve the same initial value problem. Then 

we have that the function 0 = 0 i — 02 satisfies ([!]), ([ 2 ]) and ([ 3 ]) with / = 0 
and vanishing initial data. This implies: 


if 2 < 
u+> 


< 

r mk) 2 dt 

Jo 


(47) 

< 

J o ((\\4> s 0 ) 2 + 

(II/II0) 2 )* 

(48) 


rt 1 _ 

(iioii^m 


< 

l ((l|o||^ 0 ) 2 + 

(49) 

= 

0 


(50) 


Because the norm of 


1 ^2 
ui 


) 2 = 0 that implies 0 = 0 and we have 


01 = 02 - 
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In a similar way we prove the continuity of the solution with respect initial 
data. We make the concept precise as follows. We say the solution is 
continuously stable in H X (U^, v g ) with respect initial data in H 2 (U^. v g ) 
and source functions in L 2 (U^. u g ) if for every e' > 0 there is a cq, e) 2 >c >3 
depending on ip, 7r, / such that if: 

(ll^-0lllo) 2 <*i, (51) 


7T-7r||so) <<*2 


for p, n E H 2 (Tjo, fih) and 


(II f-f\\° u+ ) 2 <s 3 


for/ E L2(U^,v g ) implies that 


II 1 i 2 < c' 

H } - 


(52) 

(53) 

(54) 


where 0 is a solution with initial data given by 0|s o = P and ^|s 0 = ^ and 
source function /. 

Choose 


<52 = 


3f i 


Then we obtain the inequalities: 


£ [j* 


(d Q p\ 2 / d°(p 


\dx c 


dx c 


Ph < 


3t i 


and 


(55) 


7r 2 — fr 2 n h ) < —- 
P “ Mi 


'S 0 


(56) 


as a direct consequence of (||<^ — ^||s 0 ) 2 < ^ and (||7r — 7r||| 0 ) 2 < 


Adding both inequalities we obtain: 


E 

a< 1 


/S 0 


d a p\ 2 (d a p 


dx c 


dx c 


Hh ) + ( / vr 2 - 7 f 2 ii h ) < ~ (57) 

'So 


2e 

3ti 


but this implies that: 
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«-£lls,) a <|( (58) 

We also have that 

(ll/-/ll“*) 2 <^ (59) 

Then applying again the energy inequality we have: 


(H0-^ni + ) 2 < 

Pu-iwhf* 

Jo 


< 

r-sj 

P m - $\\h ) 2 + (\\f 

Jo 

-/X+m 

< 

r-sj 

f h 2e e , 

J o 3t] 3ti 


< 

e 


= 

Ke 



Choosing e = we obtain the required result. 


(60) 

(61) 

(62) 

(63) 

(64) 


3 Discussion 

The equivalence between a solution of Einstein’s Field Equations and the 
background metric on which the fields propagate makes the definition of a 
singularity in General Relativity subtler than in any other physical theory. 
The most common definition of a singular-free spacetime was formulated 
by Geroch m-, a spacetime is singular-free if it is geodesically complete. 
This characterisation is well motivated as we can associate the history of 
free-falling test particles to geodesic motion. In fact, Hawking and Penrose 
[Tj showed under very general physical conditions on the topology of the 
manifold and the behaviour of the causal structure of the spacetime that 
during gravitational collapse the spacetime can not be geodesically complete. 
However, there is a gap between physical intuition and mathematical for¬ 
malism. We want to describe from a physical point of view a singularity as 
some region of spacetime that becomes more and more pathological (maybe 
by the curvature becoming unbounded) where eventually General Relativity 
is no longer adequate and some Quantum Gravity theory is needed. The 
singularity theorems make a more modest conclusion simply that spacetime 
contains a geodesic that can not be continued indefinitely. The main objec¬ 
tive of the present work is to show possible ways in which we can bridge the 
gap between intuitions and formalisms. 

In order to pursue a more complete physical description of singularities we 
proposed to define the notion of a singularity as obstruction to fields rather 
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that geodesic trajectories. Then one can talk about the field regularity of 
spacetime. In particular, in this paper the concept of classical strongly wave 
regularity was explored on spacetimes with low differentiability with certain 
timelike vector field transverse to caustics or shocks where the curvature 
is not a continuous function. The paper showed that the wave equation 
has unique solutions on the Sobolev space H l u g ) and that solutions 
depends continuously with respect initial data in H 2 (To,ph) x H 2 (Y>o,Ph)- 
In this sense the classical singularity does not disrupt the wave dynamics. 
We will further explore the field regularity of these spacetimes by extending 
this work by linking it with quantum field theory. That means to study the 
quantum wave regularity of the spacetime. Moreover, a concept of singu¬ 
lar behaviour defined in terms of quantum fields rather than geodesics will 
provide a semi-classical picture able to use concepts from Quantum Field 
Theory and General Relativity. This type of analysis can then shed light on 
the behaviour we might find in quantum gravity. 

In addition, the notion of well-posedness is related to the determinism of the 
field equations, and in the context of General Relativity this is intuitively the 
content of the Strong Cosmic Censorship Conjecture. Current approaches 
to giving a precise formulation of this difficult conjecture rely precisely on 
the degree of differentiabilty and the notion of maximal globally hyperbolic 
development HZ). It would be interesting to reformulate the conjecture 
in terms of field regularity, which can be seen as a generalisation of the 
globally hyperbolic development in the context of low differentiability. This 
will perhaps allow one to look at the conjecture in a new way. 
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4 Appendix 

In the next Proposition we show that the there is a congruence of timelike 
geodesics whose tangent vector has an essentially bounded weak derivative. 
This is a key requirement in the above discussion to make sense of the 
energy inequalities as in a general low differentiable spacetime the covariant 
derivative may be unbounded and the argument breaks down. 

Proposition 2 Let (M. g) be a Lorentzian manifold and p a point in an 
open subset P C M with compact closure such that there is a lens-shaped 
domain U satisfying the condition of Lemma [1} 

Then there exists a congruence of timelike geodesics whose tangent vector, 
T a , has an essentially bounded weak covariant derivative. 
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Summary of the proof of Lemma. 

The proof will consist of eight steps. 

1. The first step defines the geometric setting and the specific class of 
modifiers that are used. 

2. The second step focuses establishing a majorizing ordinary differential 
equation that helps to uniformly bound the norm of the tangent vectors 
of the mollified geodesics. 

3. The third step establishes a uniform time for existence of the mollified 
geodesics such that the tangent vectors of the mollified geodesics, 7 n 
are contained inside the set C (see hypothesis 5 of Lemma ([!]) ). 

4. The fourth step focuses on establishing a majorizing ordinary differen¬ 
tial equation that helps to essentially bound the geodesic connecting 
vector Y . 

5. The fifth step establishes that Y is essentially bounded. 

6. The sixth step uses the Arzela-Ascoli theorem and the bounds previ¬ 
ously obtained to show that the limit lim^^oo is well-defined and 
gives meaning to the notion of geodesic and tangent vector. 

7. The seventh step establishes the distributional nature of the weak 
covariant derivative of the tangent vector. 

8. The eighth step shows the essential boundedness of the weak covariant 
derivative. 

Remark: Using the specific class of modifiers described below in step 1 is 
important in step 4 where it allows us to establish the conjecture of Clarke. 
First Step. 

Let U be a lens-shaped domain contained in compact set U that contains p. 
Also assume there is an appropriate choice of coordinates x l = {r,x a } on 
such that 0 < r < t\ (where t\ is defined in the second step) and the 
constant values of t correspond to the spacelike surfaces S T , £ = So . (see 
figure 1) 

Now we choose a vector V E C C M 4 . Then define the vector V q as the 
vector in q G £ whose components in the vector basis of the coordinates 
{t, x"} are the same as the components of V in the canonical basis in M 4 . 
The hypothesis of the theorem only allow to define the connection as a 
distribution. In order to use the differential equations in the classical sense 
we will use convolutions and then take limits. 

The strict delta net (p n )n that we use is defined using the admissible modi¬ 
fiers described by Steinbauer and Vickers in m • We describe them below: 
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Definition 1 There exist strict delta net ( p n ) n with 

• supp(p n ) C Bi( 0) for all n £ N 

n 

• Jm 4 z 9 ” = 1 / or M n € N which are moderate, have vanishing moments 
and the negative parts have arbitrarilly small L 1 (W i ,d 4 x)-norm. This 
last condition means that 

• Vr? > 0 3Nq : f R 4 | p n | = I + 77 /or all n > Nq 

These mollifiers have the properties that if the components of g a b are in 
H} oc (M, v g ) fl Lfo c (M, u g ) then the convolution of the components g a b with 
the admissible mollifiers denoted by ( g a b)n = 9ab * Pn satisfy: 

• (9ab)n -> 9ab in H^M, Ug) n L% C (M) 

• (( 9ab)n ) 1 -» g ab in HlJM, u g ) n Tf 0C (M, u g ) for all p < 00 

• ^-'ijk((,9ab)n) ^ Tij'fc in Li oc (T/, Ug) 

• ^-fc { (Mn) ^ R) kl in £>'(«£) 

where Fy fc (( 5 a6 ) n ) is the Christoffel symbols of the metric (, g ab ) n and R) k i{{g a b)n) 
is the curvature of the metric ( g a b)n 

For convenience we are going to denote Tijk((g a b) n ) = ^tjk(n) and i n general 
all Christoffel symbols appearing will be understand as those coming from 
the convolution of the metric by the corresponding element of the delta net. 

Second Step. 

Now we define for each point q in S, a family of geodesics { 7 ™(s)}that 
satisfies: 


d 2 l q n _ T i(n)dYn j dYn k 
ds 2 ds ds 

(65) 

with the initial conditions 


si 

^0 

II 

(66) 

f <°> - * 

(67) 


We now show that there is some time si (uniform in n) such that for |s| < si 
we have ^jyf(s) £ C. We will drop the q in the notation until it is needed. 
First we define: 


l = sup{\\V\\,:VeC} 


( 68 ) 
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(69) 


r = inf{dist(-j^( 0), C c ) : q € T,} 

where dist means the euclidean distance and where C c is the complement of 
C. We first give a uniform bound 
To do this we note that: 


dr/„ 

ds 


in terms of s and the intial value. 


( 

din 

\ d 

din 


ds 

J ds 

ds 


1 d 

din 

2 ds 

ds 

d f c 

''In \ 

ds V 

ds ) 

d 

'dln\ 

ds 

v ds J 


2 

din 

ds 

din 

ds 


where u■ v is the dot product in M 4 . 
Hence 


(70) 

(71) 

(72) 


d 

din 

< 

d 

f dln\ 

ds 

ds 


ds 

l ds J 


Consider now the following inequalities: 


(73) 



( dln\ 

ds 

l ds J 


< 2sitpj = o,..,3 
= sup i= o ,.. i3 
" 32s'Up2j,/c 


d_ 

{ di l n \ 

ds 

[ds) 


-r 


i(n) dll 1 dll' 

i k ds ds 


-r 


-r 


i{n) 

jk 

i(n) 

jk 


S'UPj,k= 0,.. ,3 


din 


ds ds 


ds 


which using (73) gives: 


(74) 

(75) 

(76) 

(77) 


d_ 

ds 


din 

ds 


^^SUPij^k 


p i(n) 

1 jk 



So 


d-yn 

ds 


is bounded by the majorizing equation 


(78) 


dx 

ds 


77. 2 


\ n x 


(79) 


subject to the initial condition x'(0) < l and where X n (s) := 32 supij^ 
Then we have: 


r i(n) 

1 jk 
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£ < X n x 2 
as 

_d^ (V 

ds 


I) 5 A "« 

*( S ) ,(0) £ -/„' A » s ' 


1 


x(s) < l^oyfWw < kx{0) 

> 1 


for k > 


1—x(0) JJ A n (s')ds' 


(80) 

(81) 

(82) 

(83) 


Third Step. 

We now use the result above to show the existence of a time interval for 
which 7 remains in C. 

Notice that 


/V‘Sr’ 1 * = 

J°\ Pn * raids’ 

(84) 

= 

[ 1 [ Pn{z)T l jk (i(s') + z)dz\ds' 

J 0 J K 4 

(85) 

< 

f f \Pn(z)\\T , jk ('j{s') + z)\dzds' 

Jo j R 4 

( 86 ) 

< 

[ \Pn(z)\ f |r*fc(7(s / ) + z)\dsdz 

-'JR 4 J 0 

(87) 

< 

M(s) [ \p n (z)\dz 

J K 4 

( 88 ) 

< 

M(s)( 1 + 7 ) 

(89) 

for 77 > 0 when n > Nq where we have use (83) with k = 1 + 

so that 


Then g n (s ) = fj A n (s')ds' is a non decreasing function that satisfies lim 5 +^ 0 g n (s) 
0 for all n as a consequence of property (5) of the hypothesis of Lemma 
0 and ( |89[ ). Moreover, there is a smallest time so 7^ 0 such that k = 

1 + = i_3;(o)g”o(s 0 ) f° r some no- Then choosing s < so guarantees that 

1 + t r=ln( s i\ ds i f° r n - Notice that so can not be 0 because that 

would imply k = 1 contradicting k > 1. 

Then for s < s a integrating and taking absolute values on the geodesic 
equation (65) we have that 


ds ds 

j 6 I din 2 din k . , 

Jo 1 jk ds ds ldS 

(90) 

< 

16(r + /) 2 | Q Velds' 

(91) 

< 

16(r + l) 2 M(s)( 1 + 77 ) 

(92) 
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Now because M(s) —> 0 as s —> 0 we can make the difference as small as we 
want. So taking explicitly si < sq such that 


M{s i) < 


32(r + Z) 2 (1 + rj) 


Ensures that: 


as as 


< 2srtpi =0 ,..,3 


£'<•) - £‘<0) 

as as 


< 32(r + /) 2 M(s)( 1 + 7y) 
r 


(93) 


(94) 

(95) 

(96) 


for all s < si (where si is independent of n ). This implies that ^(s) G C. 
Then we can choos f i as sufficiently close to 0 to ensure that p is covered by 
the curves up to si 


Fourth Step. 

We now consider a 1-parameter family of initial conditions 


7n(° ,«) = ?(«) (97) 

^(0,«) = F jW (98) 


and the corresponding family of geodesics, {7n^(s)} u - Now let Y be the 
connecting vector of this family. For simplicity later on, we are going to use 
a parallel propagated co-frame on 7 $,{e a ( n ) = e^ n 'dx l } a= 0 , 1 , 2,3 coinciding 
with the coordinate basis at s = 0 so we define: 


ya ^7 V a(n) yi a(n) 

u ' du * ' u 1 


(99) 


Moreover if q(u) is a coordinate function of E i.e., q{u ) = x a , we will denote 
the connecting vector as Y a with frame components Y£. Now we have 


D_ (dY“ d \ 

Ds \ ds ds) 
d 2 Y£ 8 
ds 2 ds 

where we have use the fact that Y£ is a scalar and j is the geodesic 

equation where we have define jy := V_a_. 

ds 

Now notice that in frame components we have 




D 2 Y“ 

Ds 2 
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(102) 


D 2 YS = D (DYj M „) Yi Pef n> \ 

Ds 2 Ds\ Ds 1 a Ds ) 

- D ( dY iMn)i> , 

Ds \ ds i ds a Ds j 

D 2 Y* a{n) 8 
Ds 2 &i ds 


(103) 

(104) 


where we have use the fact that, Y£ is a scalar and n are parallel propa¬ 
gated coefficients. 

Now using the geodesic deviation equation and using (104) and (101) we 
have: 


d 2 Y a D 2 ya 

ds 2 ~ Ds 2 

(105) 

D 2 Y i i \ 

_ U l a a(n) 

Ds 2 1 

(106) 

a(n) r>i ^In’’ d^ n k b U n ) 

1 ]kl ds ds a b 

(107) 

with the initial conditions: 


4 (0) = Y& 0 )e“ (n} (0) 

(108) 

= « 

(109) 

<-0 

II 

(110) 

where the first 5 is by noting that at s = 0 we have 

- and the 


second one as a consequence of the initial alignment with the coordinates. 
Also notice that 


dJl 

ds 


=(«) 


(Vaf a j 4 n) 

\ os / 

(VaF„j eW+y Q fl Vaei") 

\ os / OS 


VaO? 

OS 

d 

V y^e ( a n) 

^4 


(in) 

( 112 ) 

(113) 

(114) 

(115) 


where (114) is the use of the torsion free condition. Now if we evaluate at 
s = 0we have the second initial condition by the following calculation 
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dYl 

ds 


(0 ) 4 ” > 


where (120) is a definition. 
Now notice that: 



(116) 

9 1 

^ V ei n) (0) fo' 0 

(117) 

„ d . 

V_9_ W- ° 
dx a ds 

(118) 

V a Vg 

dx a 

(119) 

d 

(120) 


V A e a(n) = 0 

ds 

as a consequence of being parallel propagated. 
Then in coordinates we have: 

def n) _ k(n) a^d'yj 
ds ~ ij 6k ds 
So using again a majorizing equation 


( 121 ) 


( 122 ) 


d i=cX 

ds 

with initial condition z(0) = 1, c a constant and applying similar arguments 
as the ones used in step two and three we can conclude that is uniformly 
bounded in terms of M(s i). 

It follows then that using (107) and the uniform bounds on and ^7^ 

(123) 


that we can write: 


d 2 Y“ 

ds 2 


< C<t\\Y± 


a || 

ol Moo 


where C is a suitable constant and 

cr := suplR^ 1 (i n {s))\ 
Next we consider the majorizing equation 

d 2 z 


ds 2 


= Caz 


with initial conditions 


z(0) = 1, = sup\tf\ 

ds 


(124) 


(125) 


(126) 
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where the supremum is taken with respect all indices that appear in •&. 
Notice that bounding z, implies a bound of ||h^||oo- 

Fifth Step. 

We now obtain a bound on z in terms of the curvature tensor. For that 
purpose we need the following inequality: 



•S 


< 

/ ads' 

Jo 

(127) 

= 

sup\p n * R) kl \ds' 

(128) 

= 

sup\ J p n (z)R l j k i('y(s') + z)dz\ds 

(129) 

< 

I Jo sup][Rl ^ 1 ^ s ^ + z )\ ds ' dz 

(130) 

< 

N(s ) J \p n (z)\dz 

(131) 

< 

N(s ) 

(132) 


To estimate a solution to (125) we notice that if S 2 is the first value of £ at 
which 1 2: | = 2 (possibly S2 = 00) then before s 2 


d 2 z 
ds 2 

< 


Caz 

Ca2 


(133) 

(134) 


which is a consequence of the initial conditions and the continuity of 2 . 
Now integrating two times both sides considering (127) and the initial con¬ 
ditions on z we got 


z < 1 + sup\d\s + 2 C 


N(s')ds' 


(135) 


for 0 < s < S 2 < si 

Now the right side is an increasing function that starts at zero. So there is 
a s 3 such that 


|$|s 3 + 2C p N(s')ds' < 1 
then we will have 2 < 2 up to S 3 , and hence 


W a ll < 9 

| 1 a II00 ^ 


(136) 


(137) 


in this interval. 
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Sixth Step. 

The Arzela-Ascoli theorem guarantees that given {p n } n ,a sequence of equicon- 
tinuous and uniformly bounded functions there is a sub-sequence that con¬ 
verges uniformly. 

Now the functions 




are equicontinuous and uniformly bounded. This can be seen by noting that 
the functions are defined in a bounded domain with a direct application of 
the mean value theorem in several variables and the fact that and 


are uniformly bounded. 

So by Arzela -Ascoli theorem there is a sub-sequence of { 7 ™} that gives 
meaning to the idea of a geodesic 7 with tangent vector T. 


Seventh Step. 

We now define the j component of as: 


Vi 


djn 

ds 


,<t>) = 




Vi 


drfn 

ds 


<K 


d ^d'yj 

d~/n k ( d(j) 


i( n)j d'Jn 


Ms?hr +r « * 


Ju+ ds V dx l 

£ i 


4 - r£?V) 


(138) 

(139) 

(140) 


where (140) is obtained by integrations by parts. 

Notice that the right hand converges in M as n tends to infinity for every 
4> E T>(JA ^). Hence, the expression converges in the sense of distributions to 
the distributional covariant derivative of T (see H2], p. 134). 


Eight Step. 

We now establish the essential boundedness of the weak covariant derivative 
of the tangent vector. Consider a basis for T P U^ is {^,V ( (^ T } Q= 1 , 2 , 3 . So 
any vector X can be written as a linear combination of those. 

If the i- component of X is X 1 = X a Y ( ! > + X then 
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Vx 


djn 

ds 


w d'Jn 

X xi _d_ — 
dx i ds 

Y7 | ^7 

V X“Y*-§- * + V 


ds 


djn 

x o^2?i i _a T j 
- V da dx i as 


Y~ a X7 i V°Y7 

X V y —-hA V dTjj —— 

as ^ as 

ds 

X a Vd^Y“e^ 

ds 

rW a 

X a a r M 

ds a 


(141) 

(142) 

(143) 

(144) 

(145) 

(146) 


where we have used the linearity of the covariant derivative, the torsion free 
condition and the vanishing Lie bracket between Y a , and along with the 
geodesic equation. 


Now integrating once (107) we have: 


^ = ^o)rl"Wo)) + 


“W Tji 


d'jJd'Jn 
ds 


Y*e b {n) ds' 


(147) 


Now this implies that {Xi^YY is bounded and the bound is independent 

X~~ ”” rj'Y' a 

of n. This can be seen by using (94), (122),(137) to bound and the fact 
that X a are continuous functions (any vector X is a tangent vector of a C 1 
curve) in a compact set so a bound exist. 

\i 


Now using the boundedness of [Xi^r 
the form: 


and (138) we have an estimate of 




(148) 


= //• 


where B is a co nstant and 
Moreover, (148) allows us to define (Xi^^Y as a functional over the space 
of integrable functions, L 1 (U^. fig). Now taking the limit as n 
obtain 


oo we 


,£&((*<£) (149) 

because all the bounds hold in the limit. 

So lim n _ ) . 00 (x converge in the dual space (the space of linear function¬ 
als) of integrable functions. This space is isomorphic to L°°(U^) because 
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is compact. Then G L°°(U^) and is a uniformly bounded func¬ 
tion. Finally we have that lim, woo is essentially bounded because 

each component is an essentially bounded function. 
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